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We derive the leading contribution to the four-nucleon force within the framework of chiral 
effective field theory. It is governed by the exchange of pious and the lowest-order nucleon-nucleon 
contact interaction and includes effects due to the nonlinear pion-nucleon couplings and the pion 
self interactions constrained by the chiral symmetry of QCD. The resulting four-nucleon force does 
not contain any unknown parameters and can be tested in few- and many-nucleon studies. 
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Chiral effective field theory (EFT) is a powerful tool for 
analyzing the properties of hadronic systems at low en- 
ergies in a systematic and model-independent way. It 
exploits the approximate and spontaneously broken chi- 
ral symmetry of QCD which governs low-energy hadron 
structure and dynamics. Over the past years, consid- 
erable progress has been achieved in understanding the 
structure of the nuclear force in this framework P, U, see 
for a review. In particular, two- (2NFs) and three- 
nucleon forces (3NFs) have been worked out and applied 
in few-nucleon studies upto next-to-next-to-next-to- 
leading (N^Lqi H 13 and next-to-next-to-leading or- 
ders (N^LO) in the chiral expansion, respectively. 
Parallel to these developments, significant progress has 
been achieved in the microscopic description of few- and 
many-nucleon systems. Advances in the development 
of methods, such as the Faddeev-Yakubovsky scheme. 
Green's function Monte Carlo and hyperspherical har- 
monics methods as well as the no-core shell model, 
coupled with a significant increase in computational re- 
sources allow one to obtain the binding energies of light 
nuclei with high numerical accuracy reaching 0.2% for the 
a-particle and 1-2% for heavier nuclei Three- 
and four-nucleon continuum observables are also being 
investigated. These studies open the door for precise and 
nontrivial tests of the underlying dynamics. In this con- 
text, three-nucleon forces have deserved much attention 
and were shown to be important e.g. for understanding 
the spectra of light nuclei _9] . Four-nucleon forces (4NFs) 
have not yet been explored in few-body calculations. Al- 
though the existing studies leave little room for the con- 
tribution of the possible 4NFs, see e.g. 0, this needs 
to be verified via explicit calculations. In this paper we 
make a step in this direction and derive the leading 4NF 
in the framework of chiral EFT which appears at N'^LO 
in the chiral expansion. The resulting 4NF is local and 
parameter-free. Its application in future few-body stud- 
ies will provide a challenging test of our understanding 
of the nature of nuclear forces. 

Let us first recall the structure of the effective chiral La- 



grangian for pions and nucleons, which at lowest order 
reads: 
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Here the superscripts refer to the number of derivatives 
and/or AfTr-insertions, Ft^ = 92 A MeV is the pion de- 
cay constant and gA = 1.267 is the nucleon axial-vector 
coupling. The low-energy constants (LECs) Cs and Ct 
determine the strength of the leading NN short-range in- 
teraction [US]. The Lag rangians and ^cj^rjY ^^^^ given 
in the heavy-baryon formulation with N representing a 
non-relativistic nucleon field and a denoting the Pauli 
spin matrices. Further, the SU(2) matrix U — collects 
the pion fields, C = d^" -I- ^^'u] denotes the covari- 
ant derivative of the nucleon field and = iu^d^Uu^ . 
The first terms in the expansion of the matrix [/(tt) in 
powers of the pion fields take the form, see e.g. (Ifll] : 



[/(tt) = 1- 



i 1 2 , 8q! — 1 



^ (2) 

where r denote the Pauli isospin matrices and a is an 
arbitrary constant. Notice that only the coefficients in 
front of the linear and quadratic terms in the pion field 
are fixed uniquely from the unitarity condition [/^[/ = 1 
and the proper normalization of the pion kinetic energy. 
The explicit a-dependence of the matrix U represents 
the freedom in the definition of the pion field. Clearly, 
all measurable quantities are a-independent. 

The structure of the nuclear forces can be visualized 
utilizing a diagrammatic language. Given the non- 
uniqueness of the nuclear potential, the precise mean- 
ing of diagrams is scheme-dependent. Typical examples 
include Fcynman graphs with subtracted iterative contri- 
butions, see e.g. [l(j|, irreducible time-ordered diagrams 
0, |3| or time-ordered-like graphs in the method of uni- 



FIG. 1: Disconnected diagrams at order u = 2 that lead to 
vanishing contributions to the four-nucleon force. Sohd and 
dashed hnes represent nucleons and pions, respectively. 



FIG. 2: Examples of disconnected diagrams at order i/ = 4 
that lead to vanishing contributions to the four-nucleon force. 



tary transformation jj]. The importance of a particular 
contribution to the nuclear force is determined by the 
power V of the expansion parameter Q/A. Here, Q and 
A refer to the generic low-momentum scale associated 
with external nucleon momenta or M-r; and the pertinent 
hard scale, respectively. For a A^-nucleon diagram with 
L loops, C separately connected pieces and Vi vertices of 
type z, the low-momentum dimension is given by 

V = -2 + 27V-2C + 2L + ^y,Ai (3) 

i 

with Ai — di + —rii — 2 . 

Here, is the number of nucleon field operators and 
di the number of derivatives and/or insertions of M^^. 
The interaction terms of lowest possible chiral dimen- 
sion Ai = are listed in Eq. (QJ. The power counting 
in Eq. naturally explains the dominance of the two- 
nucleon force (2NF) which starts to contribute at order 
= P, The leading 3NF formally arises at or- 
der v = 2. The corresponding contributions are known 
to vanish provided one uses an energy-independent for- 
mulation such as the method of unitary transformation 
Hi, 12, 13], see also 14, 15j for a similar conclusion based 
on the framework of time-ordered perturbation theory. 
The first non-vanishing SNFs appear at order v = 3 
At present, the 2NFs and SNFs have been worked 
out and applied in few-body studies upto orders v = A 
0,01 and 1/ = 3 0, respectively, see (3) for more details. 
Derivation of the leading 3NF corrections at order i/ = 4 
is in progress. 

According to Eq. Q , the leading 4NF is expected to arise 
at order v = 2 from disconnected tree diagrams with the 
lowest-order vertices, see Fig. ^ Using the method of 
unitary transformation, it is easy to see that both graphs 
in Fig. n yield vanishing contributions, see also [l5j for 
an earlier discussion based on time-ordered perturbation 
theory. The first non-vanishing 4NFs appear at order 
y = A where one has to take into account disconnected 
tree diagrams with one insertion of a = 2 vertex as 
well as disconnected loop and connected tree graphs with 
the lowest-order vertices depicted in Figs. HI and re- 
spectively. In addition, one needs to consider relativis- 
tic 1/m-corrections to disconnected diagrams shown in 
Fig. n We have verified that all disconnected graphs at 
this order yield vanishing contributions to the 4NF. One 



is, therefore, left with connected tree diagrams shown in 
Fig. 131 We have used the method of unitary transfor- 
mation 01 to evaluate the corresponding 4NF contribu- 
tions. To that aim we first switch from the Lagrangian in 
Eq. to the effective ttN Hamilton operator using the 
canonical formalism and then apply a suitable unitary 
transformation in Fock space in order to decouple purely 
nucleonic states from the ones which contain pions. The 
resulting Hamilton operator acts on the nucleon subspace 
of the Fock space and provides a definition of the nuclear 
forces. A rather tedious calculation along the lines of 
Ref. 01 yields the following result for the contributions 
of individual diagrams in Fig. |31 
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FIG. 3: The leading contributions to the four-nucleon force. 
Graphs resulting by the interchange of the vertex ordering 
and/or nucleon lines are not shown. 
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Here, the subscripts refer to the nucleon labels and 
Qi — Pi' ~ Pi with Pi ' and pi being the final and ini- 
tial momenta of the nucleon i. Further, qi2 — qi + 02 = 
— 03 — 54 = —(734 is the momentum transfer between the 
nucleon pairs 12 and 34. Notice that some terms in the 
sum over all possible 24 permutations of the nucleon la- 
bels are the same. Diagrams (b), (d), (g), (h), (i), (j), 
(m), (o) and (p) lead to vanishing contributions to the 
4NF. Notice further that the short-range 4NFs only de- 
pend on the LEG Ct which can be traced back to the 
fact that the Cs -vertex commutes with the other lowest- 
order vertices. Details of the derivation will be published 
elsewhere. 

The results given in Eq. |0J deserve some special com- 
ments. The easiest way to evaluate the contributions 
from graphs (e) and (f), which do not include reducible 
topologies (i. e. time-ordered diagrams with purely nu- 
cleonic intermediate states), is by calculating the corre- 
sponding Feynman diagrams. The resulting 4NFs de- 
pend individually on the parameter a which enters the 
effective Lagrangian in Eq. When they are added 
together, this a-dependence cancels out, see also [lillr^ 
for an earlier discussion of these contributions. We only 



FIG. 4: Selected contributions to the three-nucleon force. 

show explicitly the remaining, a-independent terms in 
Eq. The 37r-exchange 4NF proportional to g\ can 
be obtained by evaluating the Feynman diagram (g) in 
Fig. 13 Due to the four-momentum conservation at each 
vertex and the fact that the Weinberg-Tomozawa vertex 
contains a time derivative of the pion field, this contri- 
bution is suppressed by Q"^ jrr? and, therefore, shifted to 
higher orders. In the method of unitary transformation, 
one needs to consider not only graph (g) but also graph 
(h), where the ttttNNNN vertex of dimension Ai = 0, 
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arises in the Hamiltonian through the application of the 
canonical formalism to the Lagrangian in Eq. see 
^M- We have verified that the 4NF contributions 
from diagrams (g) and (h) in this framework cancel 
against each other, which is consistent with the result 
based on the Feynman graph technique. 

Before discussing the numerical size of the derived 4NF, 
another important comment is in order. Nuclear forces 
are not unique and can, in general, be changed via uni- 
tary transformations. We have explored a large class of 
unitary transformations which act on the nucleon sub- 
space of the Fock space and affect nuclear forces at or- 
der V — A (and higher). We found that the associated 
ambiguity in the definition of the nuclear forces is re- 
solved if one requires renormalizability of the l~loop 3NF 
contributions. Consider, for example, 37r-exchange 3NF 
corresponding to diagrams (a-c) in Fig. 0] Ultraviolet- 
divergent loop integrals resulting from 27r-exchange be- 
tween the nucleons 1 and 2 have to be renormalized by 
a redefinition of the LEG D which enters the expression 
for the 3NF at order = 3 0, , see diagram (d) in 
Fig.H 
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For this to be possible, ultraviolet-divergent pieces in the 
37r-exchange 3NF must have the same form as V^m- As 
a consequence, no terms proportional to {q^ -f AI^)^^, 

{qi + A/2)-3/2 (^g'3^+M2 + ^ri2+M2)-l with W 

being the momentum of one of the two pions exchanged 
between the nucleons 1 and 2 should appear in the 3NF. 
This restriction imposes nontrivial constraints on the 
choice of the unitary transformations and, therefore, also 
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on the form of the 4NFs corresponding to graphs (a) and 
(b) in Fig. |21 These considerations may remind one of 
the recent findings in the context of large QCD |0| . 
One should, of course, keep in mind that S-matrix ele- 
ments, which involve iterations of the potential, are not 
affected by unitary transformations. Based on the above 
arguments, we have not found any ambiguity in the form 
of the 4NFs (for the considered class of unitary transfor- 
mations). We postpone a detailed discussion of this issue 
to a forthcoming publication. 

In order to test the effects of the 4NFs in few-nucleon 
systems, explicit calculations will need to be performed. 
To get a rough idea about the size of the 4NF contribu- 
tions to e.g. the a-particlc binding energy, one can look 
at the strength of the corresponding r-space potentials. 
For V^, for example, it is given by 

(^(^M.^-80keV, (7) 

see also Here, the factor (47r)~^ arises from inte- 

grations over the relative momenta. Other 37r-exchange 
4NFs have similar strengths. The short-range terms can 
be expected to be less important due to numerically small 
values of the LEG Ct found in and (for certain cut- 
off combinations) in We remind the reader that in 
EFT without or with perturbative pions, one has Ct = 
in the limit when both NN S-wave scattering lengths go 
to infinity ^IJ. 

To summarize, we have derived the leading 4NF con- 
tributions based on chiral EFT that appear at N'^LO in 
the low-momentum expansion. Chiral symmetry of QCD 
fixes the structure and strength of the pion couplings and 
plays a crucial role in the derivation. The strength of 
the lowest-order NN contact interaction entering the ex- 
pressions for the 4NF is determined in the two-nuclcon 
system 0, so that the resulting 4NF is parameter- 
free. Applications of this 4NF to study the properties of 
the ck-particle and the spectra of other light nuclei will 
provide a challenging test of our understanding of the 
nuclear force. 
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